On Melting Temperature of Heavy Quarkonium with the AdS/CFT Implied Potential 
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The quarkonium states in a quark-gluon plasma is examined with the potential implied by 
AdS/CFT duality. Both the vanila AdS-Schwarzschild metric and the one with an infrared cut- 
off are considered. The calculated dissociation temperatures for J/tp and T are found to agree with 
the lattice results within a factor of two. 
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The heavy quarkonium dissociation is one of the important signals of the formation of QGP in RHIC. The subject 
has been explored extensively on a lattice [H, @, In the deconfinement phase of QCD, the range of the binding 
\^ ' potential between a quark and an antiquark is limited by the screening length in a hot and dense medium, which 
decreases with an increasing temperature. Beyond the dissociate temperature, X^, the range of the potential is too 
short to hold a bound state and the heavy quarkonium will melt. The lattice simulation of the quark-antiquark 
potential and the spectral density of hadronic correlators yield consistent picture of the quarkonium dissociation and 
the numerical values On the other hand, AdS/CFT duality provides a new avenue towards a qualitative or even 
CV semi-quantitative understanding of the non-perturbative aspect of a quantum field theory 0,01 . It is conjectured that 
a string theory formulated in AdS§ x S5 is dual to the N — 4 supersymmetric Yang-Mills theory (SUSY YM) on the 
boundary. In particular, the low energy limit of the classical string theory, the supergravity in AdS§ x S5 corresponds 
to the supersymmetric Yang-Mills theory at large N c and large 't Hooft coupling A = N c gY M . Remarkable success 
has been made in the application of the AdS/CFT duality to the physics of quark-gluon plasma (QGP) created by 
RHIC, even though the underlying dynamics of QCD is very different from that of a supersymmetric Yang-Mills 
theory [1, 0, It is the object of this brief report to calculate Td using the heavy quark potential of N = 4 SUSY 
YM extracted from its gravity dual d, [l(| H3 • 

Although the potential model applies only in the non-relativistic limit which is not the case when the 't Hooft 
coupling, A, becomes too strong, it can be justified within the lower side of the range of A used in the literature to 
compare AdS/CFT with the RHIC phenomenology, i.e. 5. 5 < A < Qir. The upper edge is obtained by using N c = 3 
and the QCD value of gyivi at RHIC energy scale (gyu/i^) — 1/2) and the lower edg e is based on a comparison 
, between the heavy quark potential from lattice simulation with that from AdS/CFT 

We model the quarkonium, J/ip or T, as a non-relativistic bound state of a heavy quark and its antiparticle . The 
. £h ' wave function for their relative motion satisfies the Schrodinger equation 

- 7 L V ^ + VcS ( r )yj = -Evh (1) 

where m = M/2 is the reduced mass with M the mass of the heavy quark and E(> 0) is the binding energy of the 
bound state. Because of the screening of QGP, the effective potential energy has a finite range and is temperature 
dependent. The dissociation temperature of a particular state is the temperature at which its energy, —E, is elevated 
to zero. 

The free energy of a static pair of qq separated by a distance r at temperature T is given 

e -iF(r,T) = tr < WHL + )W(L_) > 

tr< Wi(L+) >< W(L_) > [ ' 

where L± stands for the Wilson line running in Euclidean time direction at spatial coordinates (0,0, ±5?") and is 

T 1 

closed by the periodicity. We have W(L±) = Pe~ l f° dt - 4 o(t,o,o,± 5 r) w ^h Aq the temporal component of the gauge 
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potential subject to the periodic boundary condition, Ao(t + ^ , r) — Ao(t,r). The trace here is over the color indexes 
and < ... > denotes the thermal average. The symbol P enforces the path ordering. The corresponding internal 

energy reads U(r, T) == —T 2 -§^ ( F ^ T ^ ■ Two ansatz of the potential model have been explored in the literature [1]: 

the _F-ansatz which identifies V e g. of |T]) with F(r, T) and the JJ-ansatz which identifies V s. with U(r, T). The lattice 
QCD simulation reveals that the U ansatz produces a deeper potential well and thereby higher because the entropy 
contribution is subtracted. This remains the case with holographic potential as we shall see. 

According to the holographic principle, the thermal average of a Wilson loop operator W(C) = Pe~ l §c in 
4D N — 4 SUSY YM at large N c and large 't Hooft coupling corresponds to the minimum area S m i n . [C] of the string 
world sheet in the 5D AdS-Schwarzschild metric with a Euclidean signature , 

ds 2 = n 2 T 2 y 2 (fdt 2 + dx 2 ) + -^jdy 2 , (3) 

y j 

bounded by the loop C at the boundary, y — > oo, where / = 1 — Specifically, we have 

tr<W(C)>=e- VXs ^ [C] . (4) 

For the numerator of {2}, C consists of two parallel temporal lines (t,0, 0, ±£) and the string world sheet can be 
parameterized by t and y with the ansatz x\ = X2 = and 23 = a function of y. The induced world sheet metric 
reads 



ds 2 = Tr 2 T 2 y 2 fdt 2 



2 T 2 2f dx 3\ 2 1 

7T T y — 



V dy J ir 2 T 2 y 2 f 
Minimizing the world sheet area (the Nambu-Goto action) 



dy 2 . (5) 



S[C] = (ttT)JJ dt^ d2/yi + ^TV/(^) 2 (6) 
generates two types of solutions @ [l(| [U- One corresponds to a single world-sheet with a nontrivial profile 23(2/), 

x 3 = ±nTq f ^ _= (7) 

Jyc Viv i)(y -J/c) 



where q is a constant of integration determined by the boundary condition 

2? r 00 rfy 



(8) 



with = 1 + q 2 . The corresponding value of vA5[C] is denoted by 7i. The other solution consists of two parallel 
world sheets with £3 = ±| extending to the black hole horizon and the corresponding value of y/\S[C] is denoted by 
1% . The latter solution corresponds to two non-interacting static quarks in the medium and is equal to the denominator 
of @ 

The free energy we are interested in reads 

F(r, T) = Tmin(1, 0) (9) 

where 



i = h-h = V\ 




1 - Vc 



(10) 



Inverting eq. (|HJ) to express q in terms of r and substituting the result to (fT0|) , it was found that the function I consists 
of two branches, The upper branch is always positive and is therefore unstable. The lower branch starts from being 
negative for r < ro and becomes positive for r > tq. Both branches joins at r = r c > ro beyond which the nontrivial 
solution ceases to exist. Numerically, we have tq ~ and r c ~ Introducing a dimensionless radial coordinate 

p = irTr, we find that 



F(r,T) = -^<p(p)0( Po -p), (11) 
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where a — ^f^rjV'A ~ 0.2285\/A, and <j)(p) = —pi /(net) is the screening factor . We have (f>(0) = 1 and 4>(po) = 

with po = 0.7541. 

The small p expansion of <f>(p) is given by 

r 4 (I) 3r 8 (±) 

M-i-^+y£f+o<f). (i2) 

On writing the wave function tp(r) = ui(p)Yi m (f), the radial Schrodinger equation for a zero energy bound state reads 

d 2 ui , 2dui [1(1 + 1) 



dp 2 p dp 



P 2 



V 



ui = (13) 



with V = MV cff ./(n 2 T 2 ). We have 



for the F ansatz and 



V = -^-cf>(p)9(p - p) (14) 
PoP 



PoP 



<I>(P)-P{f)\0{po-P) (15) 



for the U ansatz, where r? = \ ap '^ f 



Note that the potential of the U-ansatz jumps to zero from below at p = po, since the derivative of <j)(p) is nonzero 
there. For both ansatz, and the case with an infrared cutoff discussed below, the solution to (fl3|) is given by 

ui = const. p - ' -1 (16) 

at p > po and by 

ui = const. p l (17) 

near the origin. The threshold value of r\ at the dissociation temperature, r/ ( i, is determined by the matching condition 
at p = po, 

|(^)l P=Po -=0. (18) 
It follows from 77 that the dissociation temperature is given by 

T* = ^ = -^VXM. (19) 

r4 (3) % 

It is interesting to observe that the extrapolation of the first two terms of (|12l) vanishes at p = p' = 47r 3 /r 4 (|) — 
0.7178, which is very close to the exact zero point, and the third term of (fT2|) remains small there. This suggests that 
the screening factor (j)(p) can be well approximated by a linear function 

4>(p) ~ 1 - 4- (20) 
Po 

with po = ^(po + Po) — 0.7359. The effective potential V e g, is then approximated by a truncated Coulomb potential. 
We have 

v = -^(i-P-)e{ Po -p) (21) 

Pop V Po / 

for the F-ansatz and 

V = -^0(po-p) (22) 
Pop 



for the U-ansatz, where the over bar of po has been suppressed. 
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ansatz 


J/^(1S) 


J/i>{2S) 


j/i>{ip) 


T(1S) 


T(2S) 


T(1P) 


F 


67-124 


15-28 


13-25 


197-364 


44-81 


40-73 


U 


143-265 


27-50 


31-58 


421-780 


80-148 


92-171 



TABLE I: Td in MeV's under the truncated Coulomb approximation. The lower value of each entry corresponds to A = 5.5 
and the upper one to A = 6tt. 



The radial wave function of the F-ansatz under the truncated Coulomb approximation can be expressed in terms 
of the confluent hypergeometric function of the 1st kind for p < po, i.e. 

ui =p\F 1 {l + l- 7 J-;2l + 2;2r)-^). (23) 



Po 



The matching condition (|18[) yields the secular equation for rj, 



v \ 1 F 1 (l + 2-%;2l + 3;2r ) ) 



As r\ is reduced from above, we expect the bound states of the same I to melt successively. Therefore the hrst positive 
root corresponds to the minimum binding strength for a bound state of angular momentum I and the 2nd one to 
the threshold of the first radial excitation. Knowing the values of these 77's, the disassociation temperature can be 
computed from the formula (|19[) . For example, the threshold 77 of the IS state, r/is — 1.76, which implies that 

T d ~ 0.0173VAM. (25) 

In case of the (7-ansatz under the same approximation, we find that 

u l = ^-J 2l+1 (2 V f^) (26) 

VP v V po ' 

for p < po with J u (x) the Bessel function. The secular equation for r\ reads 

2* + l-,^±«=0. (27) 

We have i]is = 1.20 and 

T d ~ 0.0370\/AM. (28) 

Numerical results for the dissociation temperature of quarkonium are tabulated in table 1, where we have used the 
mass values M = 1.65GeV, 4.85 GeV for c and b quarks . The errors caused by the truncated Coulomb approximation 
are within 4 percent, as is shown by the numerical solution to the Schrodinger equation of the exact potential. 

Because of the conformal invariance at quantum level, there is no color confinement in N = 4 SUSY YM even 
at zero temperature. In order to simulate the confined phase of QCD at low temperature, an infrared cutoff has to 
be introduced that suppress the contribution of the AdS horizon. Two scenarios explored in the literature are the 
hard- wall model and the soft- wall models 0, [H| ■ The gravity dual of the de-confinement transition is modeled as 
the Hawking-Page transition from a metric without a black hole at T < T c to that with a black hole at T > T c . Heavy 
quark potential and the meson dissociation temperatures calculated with the hard-wall model arc identical to what 
calculated above with the vanila AdS-Schwarzschild metric. 

In case of the simplest soft- wall model (|15j). a dilaton is introduced and the gravity dual of the free energy is the 
given by 



F = - T m^J dX ] P0 d ^V9(R~12), (29) 

where c is determined by the p-mass and the transition temperature is predicted as T c ~ 0.2459m p (l6| . A variant of 
the soft-wall scenario proposed in ref . [13, EB| , admits a string frame metric with a conformal factor, i. e. 

bz 2 

ds 2 = —(.fdt 2 + dx 2 + f- l dz 2 ), (30) 
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ansatz 


J/i> 


T 


F 


NA 


235-385 


U 


219-322 


459-780 



TABLE II: Td in MeV's for the IS state with the deformed metric. "NA" means that there is no bound states above T c and 
the entry for the T with U ansatz and a — 6n is taken from the table I, since no significant increment is observed. 



ansatz 


J/ip (holographic) 


J/i/>(lattice) 


T (holographic) 


T (lattice) 


F 


NA 


1.1 


1.3-2.1 


2.3 


U 


1.2-1.7 


2.0 


2.5-4.2 


4.5 



TABLE III: The ratio Td/T c for the IS from the holographic potential and that from the lattice QCD 



The value of b = 0.184GeV 2 was obtained by fitting the lattice simulated transition temperature T c = 186MeV [HI]. 
Following the steps from ([3]) to (fTTj)) ,we can calculate the effective potentials at both ansatz [2 lj. To determine the 
dissociation temperature in this case, we have to solve the Schrodinger equation (|13|) numerically with the numerically 
calculated heavy quark potential for both ansatz, since the truncated Coulomb potential no longer approximates well. 
The modified dissociation temperatures are tabulated in the table 2, which show an significant increment in the 
vicinity of T c . The comparison between the ratios Td/T c we calculated here with that obtained from the lattice QCD 
is shown in table 3[2l|. 

In summary, we have calculated the dissociation temperatures of heavy quarkonia using the static potential implied 
by the holographic principle with both the vanila AdS-Schwarzschild metric and the one with an infrared cutoff. While 
estimations of Td have been made in the literature based on various holographic models [l9l I20I [22| , a determination 
of Td from the Schrodinger equation within the same framework remains lacking. Our work is to fill this gap. The 
authors of [l|| gave an order of magnitude estimate of the dissociation temperature relying on the screening length 
only. The author of [2(| generalized the spectral analysis of the light mesons to heavy mesons. Their criterion for the 
dissociation, however, appears slightly ad hoc and is again independent of the coupling. Both the screening length and 
the coupling strength ought to affect the heavy quaronium binding. Carrying out the analysis of the potential model 
inspired by the holographic principle to the same extent of that of QCD will address both contributions, especially 
the consistency of the range of the coupling constant extracted from the jet quenching with the heavy quarkonium 
physics. Also a detailed bound state calculation enables us to assess the validity of the non-relativistic approximation 
behind the potential model. On comparing our results with that from the lattice simulation we found that our 
ratios Td/T c extracted from the modified AdS-Schwarzschild metric ([50]) are lower than the lattice ones within a factor 
of two. That the increment in Td/T c from the F-ansatz to the U-ansatz is about a factor of two is similar to what 
reported in One has to bear in mind that the lattice results reviewed in Q were extracted from a pure SU(3) 
gauge theory without a matter field. On the other hand the matter field contents of N — 4 SUSY YM are larger than 
that of QCD with light quarks. It is possible that the additional screening effect of the matter field in N = 4 SUSY 
YM makes the heavy quarkonia more vulnerable and thereby lowers the dissociation temperature. This is consistent 
with the observation that the potential well becomes wider in the metric with the IR cutoff introduced in [It], HI] 
since some of degrees of freedom becomes massive. 
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